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We compute the two-particle matrix element and Josephson tunneling amplitude in a two- 
dimensional model of topological superconductivity which captures the physics of the doped Mott 
insulator. The hydrodynamics of topological electronic liquid consists of the compressible charge 
sector and the incompressible chiral topological spin liquid. We show that ground states differing 
by an odd number of particles are orthogonal and insertion of two extra electrons is followed by the 
emission of soft modes of the transversal spin current. The orthogonality catastrophe makes the 
physics of superconductivity drastically different from the BCS-theory but similar to the physics 
of one- dimensional electronic liquids. The wave function of a pair is dressed by soft modes. As 
a result the two particle matrix element forms a complex d-wave representation (i.e. changes sign 
under 90° degree rotation), although the gap in the electronic spectrum has no nodes. In contrast 
to the BCS-theory the tunneling amplitude has an asymmetric broad peak (much bigger than the 
gap) around the Fermi surface. We develop an operator algebra, that allows one to compute other 
correlation functions. 



PACS number(s): 74.20.Mn 
I. INTRODUCTION 

The phenomenon of "superconductivity" — the exis- 
tence of a metastable quantum state with a current in 
a macroscopical system — manifests itself as a particular 
set of correlations in the ground state: 

(i) Meissner effect 
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here jj_ is the transversal current V • jj_ = 0, and A = 

/ 2 \l/2 

I e i J is the London penetration depth and p s is 

the superfluid density; 

(ii) Gap in the electronic spectrum: singularity of the 
one-particle Green function u = fi(p) closest to the ori- 
gin in the to plane 

G(w,p) = <4(w,p)c lf (w J p)). (2) 

(iii) Anomalous expectation value 

(N + 2|c t T (r 1 )cj(r 2 )|A0 = A(n - r 2 ) ? 0. (3) 

The matrix element between the ground states of the 
system with N and N + 2 particles does not vanish in a 
macroscopic system. It gives rise to Josephson tunneling, 
but not necessarily to an off-diagonal long range order. 

These three correlations describe very different sides of 
the phenomenon: (i) hydrodynamics of an ideal liquid, 
(ii) a gap for one-particle excitations, and (iii) a two- 
particle matrix element. Nevertheless, due to the mean 
field character of the BCS theory all of them turn out 
to be essentially the same — all these quantities can be 



expressed in terms of one complex function A(ri — r 2 ). 
This misleading "advantage" of the BCS theory often al- 
lows one to draw conclusions about the gap function (ii) 
by looking at the matrix element (iii) and vice versa. 
However, the gap, Josephson current, and the penetra- 
tion depth are essentially different quantities: the first 
one characterizes the spectrum, the second is a matrix 
element, determined also by the phase of the wave func- 
tion, while the third measures transversal current-current 
correlations. 

In an electronic liquid where the interaction is strong, 
one also expects to see a difference between dissimilar im- 
plementations of superconductivity This differ- 
ence becomes dramatic in the topological (anyon) mech- 
anism of superconductivity, where the entire effect of su- 
perconductivity is due to peculiar quantum phases of 
wave functions of the ground state and low energy ex- 
citations. 

In this paper we study tunneling in two-dimensional 
topological superconductors and show that: 

(i) like in the BCS-theory, although due to entirely 
different physics, the Josephson tunneling amplitude is 
proportional to the equal time matrix element (||), and 

(ii) the phase of the pair wave function (||) depends on 
the direction of the vector ri — r 2 and forms the complex 
irreducible d-wave representation of the group of rota- 
tions of the plane. The phase difference between order 
parameters in points (1,2) and (3,4) is the twice the angle 
ip between vectors ri r 2 and r3 — r 4 
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At the same time the gap function Q(p) = yAo + v 2 p 2 



I 



has no nodes. 

Moreover, the pair wave function can also be written as 
a product of the BCS wave function and matrix element 
of the soft modes of the transversal spin current in the 
coordinate space, averaged over the Fermi surface with 
the weight e iarg ( k ) 



A(r) ~ e -«»s« Idk f D kf (r)A BC s(kf,r), 
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A B cs(k/,r) 



kdk e 



ifck f 



((e(k)-^ + Ag)Va 



(5) 



(0) 



is the BCS wave function of two particles with a relative 
momentum directed along = k//fc/ and the propaga- 
tor of soft modes is given by 
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The factor 2 in the angular dependence reflects the 
double degeneracy of zero modes and eventually the spin 
of the electron. The complexity of the tunneling am- 
plitude is a result of the violation of the time-reversal 
symmetry — an inherent feature of layered topological flu- 
ids with an odd number of layers. [In more realistic, 
anisotropic three-dimensional systems or in the layered 
system with even number of layers the time reversal sym- 
metry is restored due to alternating signs of parity break- 
ing in consecutive layersti-slii] 

Below we consider the incommensurate case, where the 
dispersion e(k) — fx ~ w(]kj — kf) is generic and the Fermi 
surface is a circle (Sec. Ill), and separately the commen- 
surate case on the square lattice near half filling (Sec. 
0). 

Results for the commensurate case have been reported 
in Ref.El. In this case the Fermi surface consists of 
four pockets around k/ = (±7r/2, ±7r/2) with dispersion 
e(k)-n ~ Efc /= (±f ,±f) v (\ k \~ k f)- Formulae (|J|@) in 
this case give 



X + iY 
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- 4arg(kf)+lkfR A BC s(R), (9) 



k f =(±f,±f) 



where (X, Y) are integer coordinates of the lattice vector 
R = r — r'. The structure of the the pair wave function 
in momentum space is shown in the Fig.[TJ. Although 
the angular dependence in (||) is more complicated than 
in the incommensurate case (ra), it also realizes the irre- 
ducible d-wave representation of the group of rotations 
of the square lattice. In particular it changes sign under 
90° rotation A(-Y,X) = -A(X,Y). 

The angular dependence of tunneling kipanyon super- 
conductors has been also studied in Refs.ujj. 




FIG. 1. The structure of the pair wave function in mo- 
mentum space: A(k) is maximal around four Fermi points 
k/ = (+^,+^) and consists of four similarly oriented unit 
vortices around each Fermi point. The phase of each vortex 
is relative to the direction of the k/ as in the eq. (jif]). It 
changes sign under 90° rotation A(— k v , k x ) = —A(k x ,k y ). 

One of the motivations q£ this work was the corner- 
SQUID-junction experimcnto in which the relative phase 
of tunneling amplitudes (|J) on different faces of a single 
crystal of YBCO has been measured. It is found to be ir 
in accord with d-wave superconductivity and apparently 
in agreement with the topological mechanism (|). 

These and many other anomalies in tunneling, trans- 
port, and the photoemission spectrum are due to a sin- 
gle phenomenon: the orthogonality catastropheO. In a 
strongly interacting environment the addition of an ex- 
tra particle to the system drastically changes its ground 
state. As a result the overlap between the ground states 
differing by an odd number of particles vanishes in the 
macroscopical system. Contrary ground states which dif- 
fer by an even number of particles (|^) are not orthogonal 
(compare to Rem). 

Tunneling involves a correlation between particles with 
different spins. Another goal of this paper is to develop 
a method to treat the spin correlation. 

This paper consists of four distinct parts. In Sec. || 
we discuss the Josephson tunneling in the presence of or- 
thogonality catastrophe. Then, in Sec. [II we give a phe- 
nomenological description of the topological mechanism 
of superconductivity. We describe hydrodynamics in Sec. 
[II A and include Fermi surface in the consideration in 



Sec. [II D . In Sec. IIIB we present a field theory model 

§) which yields the phcnomcnological picture of Sec. 
We calculate tunneling amplitude and off-diagonal 
matrix elements for the incommensurate case in Sees. 



1VB,[VC a nd fo r the nearly half-filled (commensurate) 
case in Sec. VC Finally, in the Appendix we describe a 
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way to derive the model of topological superconductivity 
from a t-J model, a canonical model of correlated elec- 
tronic systems. Although this section is of interest in its 
own right, we present it separately as an appendix for 
the sake of continuity of the text. 



II. TUNNELING AND ORTHOGONALITY 
CATASTROPHE 

Josephson tunneling is peculiar in the presence of 
an orthogonality catastrophe. The Josephson current 
through the junction (at zero bias voltage) between two 
superconductors, of which-one is conventional, is given 
by the well known formulaQ 
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where T^ p is a transmission amplitude of the junction, 
Ek = Ai | 2 + t\ is the spectrum of the conventional 
superconductor and F (p, id) is the spectral function of 
the superconductor of interest 

x{N+l\c <T ,(-p)\N + 2)6{u-e p ). (11) 

Here the sum goes over all quantum states with one ex- 
tra particle (e p is the energy of an intermediate state). 
If the spectrum is symmetric with respect to adding or 
removing a particle, i.e., F(p, id) is an odd function of ui, 
we obtain 



/=-2In^T fcp T_ fc _ p ^ 
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In the following we assume for simplicity that the trans- 
mission amplitude Tk p is strongly peaked at k, p close to 
the direction normal to the junction. This simplification 
should not change the phase dependence of the Josephson 
current although it can change the value of the critical 
current. Assuming that the gap in the superconductor 
of interest Ao is bigger than the one in the conventional 
superconductor Aq > Ai we obtain: 



JHTlVolAxI^ / 
i_ i Jo 



did 



■sin(0 o - <p(p n ,u))\F(p n ,u)\, 



(13) 



where p„ is the component of momentum normal to the 
surface of the junction and averaging over p„ is deter- 
mined by the actual form of the transition amplitude. 
In this formula (f>o and vq are the phase and density of 
states of the conventional superconductor and </>(p, uS) is 



the phase of the F-function ([LI]) of the superconductor 
of interest. 

In the BCS theory the F-function has a peak at the 
gap oj ~ Ao, such that the width of the peak is also of 
the order of the gap: 



dpF B cs(p,u) 
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(14) 



The peak gives the major contribution to the integral 
( |l3| ) . It selects a characteristic energy of the intermediate 
state e q ~ Ao and gives rise to the traditional BCS pic- 
ture of tunneling: a pair decays into two electrons while 
tunneling, so electrons tunnel independently. Short time 
processes (at id ~ e/) do not contribute to the integral 

The situation is drasticallyjdifferent in the orthogonal- 
ity catastrophe environment!!! We will show that in a 
topological superconductor an individual matrix element 
(N\c a (p)\N + 1) acquires an additional factor 1/Lkf, 
where L is the size of the system and therefore vanishes 
in a macroscopical sample, i.e., the ground states with TV 
and N + 1 particles are almost orthogonal. Nevertheless, 
the tunneling, i.e., a matrix clement between states with 
N and N + 2 particles is nonzero due to a large num- 
ber of low-energy intermediate states contributing to the 
sum ([ll]). A result of this is that the spectral function 
( |IT| ) acquires an additional factor td/ej. In contrast to 
BCS at oj ^> Ao we have 



dpF M ~^(^)~^. 



(15) 



Therefore, the characteristic scale of spectral function is 
shifted to the ultraviolet and becomes of the order of the 
Fermi energy — much larger than the scale of the gap [A 
simi lar ph enomenon occurs in momentum space. In the 
Sec. |IV C we show that the pair wave function A(/c) has 
a long tail well away from the Fermi surface (see also Fig. 
H)]: the integral (|l3| ) is saturated by ui ~ e/ ^> Ao. This 
means that a pair remains intact during the tunneling 
and the tunneling amplitude is determined by the equal 
time value of F(p,t = 0), i.e., by the matrix element 
A(r) of an instantaneous creation of a pair (||). Let us 
notice that the correction to the spectral function in the 
cq. ( |l5| ) — (id/ef) a at a = 1 is just marginal. Were a 
be less than 1 the time of the tunneling would be of the 
order of Aq 1 (A similar phenomenon has been discussed 
by Chakravarty and Anderson in the context of interlayer 
tunneling in cuprate superconductorscl) . 

In the coordinate representation the Josephson current 



I - Ime^°A(fc/n) ~ sin(0 Q - 0(fc/n,w = 0)). (16) 



In the corner-SQUID-junction geometryoQ one can di- 
rectly measure the difference of phases of A between two 
faces of the superconducting crystal. 
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The instantaneous character of tunneling and the 
power laws (|l^) are known from one-dimensional elec- 
tronic systems where the orthogonality catastrophe 
comes to its own. In addition, in 2D it also leads to the 
angular dependence of the two-particle amplitude (||J5|). 

The direct tunneling current is also strongly affected 
by the orthogonality catastrophe: 

I^ = 2\T\^f g/^I-^k). (17) 

Elsewhere we will show that due to the orthogonality 
catastrophe the one-particle Green function acquires a 
branch cut, rather than a pole on the threshold. As a re- 
sult the density of states in a topological superconductor 
just above the gap is suppressed by the factor \Jlo — Ao. 
This leads to a suppression of direct current close to the 
threshold bias voltage. With a logarithmic accuracy we 
will obtain almost linear I-V behavior 

J dlr ~(y-A ). (18) 

This is to be compared to the BCS-theory result: idir ~ 
(V-A )^. 



III. TOPOLOGICAL MECHANISM OF 
SUPERCONDUCTIVITY 

A. Hydrodynamics 

The mechanism of topological superconductivity is 
a generalization to higher dimensions of the Peierls- 
Frohlich phenomenon, ,-kpown from one-dimensional 
electron-phonon systemsotJ. It inherits major features 
from electronic physics in one dimension. We attempt to 
present it here without addressing a particular model but 
accepting some minimal assumptions. Later in the Ap- 
pendix we develop a microscopical model to justify and 
illustrate the phcnomcnological picture. 

(i) Zero Modes and Topological Instability. Let us con- 
sider an electronic liquid where the interaction between 
electrons is mediated by an electrically neutral bosonic 
field, that can form a point-like spatial topological con- 
figuration (soliton). Let us suppose, that in a sector with 
zero topological charge the electronic spectrum has a gap 
to. Assume now that in the presence of a static soliton 
the electronic spectrum differs from an unperturbed one 
by an additional state just at the top of the valence band 
or within the gap — a so-called zero mode or a midgap 
stateEj. If the zero mode is separated from the spec- 
trum, its wave function is localized around the core of 
the topological defect. In case when the level is attached 
to a band, the wave function decays as a power law away, 
from the center of the soliton. A general argumentts 
suggests that the midgap state always has an even de- 
generacy. This degeneracy leads eventually to a proper 
flux quantization and below it is assumed to be twofold. 



Now let us add an even number of extra electrons with 
a concentration 5 into the system. They may occupy a 
new state at the Fermi level of the conduction band. It 
costs the energy of the gap plus the Fermi energy m + /iy 
per particle, where pf is a chemical potential. Alterna- 
tively, the system may create a topological configuration 
and a number of zero modes in order to accommodate 
all extra particles. The energy of this state is the soli- 
ton mass plus exponentially small corrections due to the 
interactions between zero modes. If the latter energy is 
less than p f then every two extra electrons added to the 
system create a soliton and then completely fill a zero 
mode, rather than occupy the Fermi level of the state 
with zero topological charge. As a result the total num- 
ber of solitons in the ground state is equal to half of the 
total number of electrons in the system. 

Formally it means that, contrary to the Landau Fermi- 
liquid picture, the expansion of the energy in small 
smooth variation of chemical potential Sp(r) has a non 
vanishing linear term in Sp: 

5E{p) = J Sp{r)pdr + f 6p,(r)K(n,T - r')^-^ drdr' 

J 27T 

+ 0(5p 2 ). (19) 

Here p is the electronic density, F(r)/2ir is a density of 
a topological charge and K(p,r) is some kernel. The lin- 
ear term in chemical potential is known as Chern-Simons 
term. 

The minimum of energy is achieved if the variation 
of density is followed by the variation of the topological 
charge 

p{v)-p = J K(»,r-r>)^dr>. (20) 

While doping, electrons create and occupy zero mode 
states to minimize their energy, thus giving a non zero 
value to the topological charge. Due to the twofold de- 
generacy of zero mode states 

J K{p,v)dv = 2 (21) 

so a number of extra particles S gives rise to a flux 

5 = 2 [ ^-dr. (22) 

J 27T 

We refer to this phenomenon as topological instability. 
Once zero mode states are occupied, residual interaction 
between them lifts the degeneracy, so that zero mode 
states form a narrow band. This band is always com- 
pletely filled and is detached from the rest of the spec- 
trum by some gap Ao- 

This is already sufficient to conclude about supercon- 
ductivity — the chemical potential always lies in a gap. 
The following arguments are borrowed from Frohlich's 
papertl. The position of the topological excitation is not 
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fixed relative to the crystal lattice. Therefore, a pair 
of electrons bound to a topological excitation can eas- 
ily slide through the system (and therefore carry elec- 
tric current). It slides unattenuatedly, since the state is 
completely filled and is separated by the gap from the 
unoccupied electronic states. As a result, the low en- 
ergy physics of density fluctuations is described by the 
hydrodynamics of a liquid of zero modes: 



2p 



(23) 



where /?, vq are the average density and the velocity of the 
sound mode, j = pv is the current and v is the velocity 



[p(r),v(r')] = -iV r J(r-r'). 



(24) 



In dimensions higher than one, where the pinning effects 
are not that important, the eqs. ( p3| , p4| ) already imply 
the Meissner effect (0) and superconductivitylijo. 

It is convenient to re-parameterize densities and cur- 
rents in terms of charge displacement u: 



P~ P 



Vu, 



dtM, 



so that 



[ji(r), u/c(r')] = ~ipkkK T ~ r ')- 



(25) 



(26) 



Let us stress that this mechanism is very different from 
the mechanism where an electronic pair is localized by a 
polaron. In contrast, in the topological mechanism the 
electric charge of a pair is only partially localized at the 
core of the topological soliton. Although the number of 
zero modes is equal to the topological charge of the soli- 
ton, a part of the charge is smoothly distributed through- 
out the rest of the system. 

(ii) Single-valuedness and gauge invariance. The elec- 
tronic wave function is single- valued. This important 
property may be reformulated in terms of gauge invari- 
ance. In two spatial dimensions a topological density can 
be described by means of a gauge field A: 



F(r) = V x A. 



(27) 



The single-valuedness of the electronic wave function 
means that electrons are neutral with respect to this 
field: the electronic operator does not transform under 
the gauge transformation A — > A + VA, i.e., 



[A i (r),u k (r')]=0; [A l (r), Jk (r')}=0. 



(28) 



Let us note that, despite of what eq. (|20| ) suggests, the 
vector potential Ai is not equal to the displacement e tk Uk 
but differs from it by a singular gauge transformation. 

In fact the wave functions of zero modes represent con- 
formal blocks of some conformal field theory. 

(Hi) Incompressible chiral spin liquid. The hydrody- 
namics of charge degrees of freedom alone is not suf- 
ficient to draw conclusions about fermionic matrix ele- 
ments. They are determined also by the distribution of 



electronic spin within the zero mode. Unless there are 
solitons which create zero modes in spin sector, the com- 
pressibility with respect to density modulations generally 
results in incompressibility with respect to spin modula- 
tions. Moreover, the spin liquid in the singlet sector is 
a topological liquid. In the next section we show that 
on a very general basis, spin displacements u s defined as 
S 3 = — 1/2 Vu s obey anomalous commutation relations 



Mr) x u s (r')] = 2i7T<5(r-r') 



(29) 



in addition to standard [?'»,« (r),Mjfe )a (r')] = —ipSikS(r — 
r'). This property determines rotational properties of a 
local singlet it has I = 2 angular momentum — and even- 
tually determines the phase of the tunneling amplitude 
(§• 

By combining charge and spin parts of the hydrody- 
namics (^3[ ^) we obtain the hydrodynamics of topolog- 
ical superconductivity: 



C = C C + C S1 

C c = U{d t uf 
2p 



«o(Vu) 2 ), 
1 



C s = 2iru s x d t u s + -L((9 t u s ) 2 - d 2 (Vu s ) 2 ) 
2p 



(30) 
(31) 

(32) 



The hydrodynamics consists of two independent fluids: a 
compressible charged liquid 



(u"(w,k),u"(-w, -k)) = 
(u-V.kJ.u-^-w.-k)) = 



p 



p 



(wok) 



(33) 



and an incompressible topological (chiral) spin liquid 



(«!(w,k),!il(- u .- k )) = 



<^(c,k),^(-c,-k)) = -^(1 - ^ 
(uj!(£j,k),uj-(-w,-k)) = — 



(34) 



where $ 



4-7T and 



u", u 



are the longitudinal and 



ki 



a" 



transversal parts of the displacement Uj(k) 

e fc fcj u^. Eqs. p3] ) imply compressibility of the charge 
fluid and the Meissner effect. The absence of the sound 
mode in the spin liquid (Q) is the consequence of the 
topological mass generated by the Chern-Simons term in 

(iv) Edge Spin Current. One of the direct consequences 
of the chiral nature of the spin liquid is that the spin liq- 
uid generates spin edge current. Indeed, the spin sector 
of the hydrodynamics ([52]) is equivalent, to the hydrody- 
namics of the FQHE fluid (see e.g. JilMLS). Similar to the 
FQHE the spin excitations are suppressed in the bulk 
but develop a spontaneous spin edge current with the 
level k = 2 current algebra. Let us stress that the spin 
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edge current is the only hydrodynamical manifestation of 
spontaneous parity breaking. Contrary to a number of 
claims scattered through the literature, the spontaneous 
parity breaking is invisible in the charge sector even for 
the systems with an odd number of layers. 



B. The field theory 



Let us realize the commutation relation 
by means of an electron operator: 



§E8 29) 



p = c f c; 



S 3 = - C Vc; j 



j = —ic)Vc\ 
1 . 
2 



ic<<j 3 Vc. 



(35) 



The topological constraint ( j22| ) reduces the Hilbert space 
of the low-energy sector of the theory. Instead of using a 
projected electronic operator c a we introduce a "spinned" 
fermion operator ?/v in which terms the field theory which 
takes topological constraint ( p2| ) into account is local. We 
discuss the relation between the physical electron opera- 
tor c a and the unphysical gauge non-invariant operator 
ip a in the Sees. [III D||VA. 



In tfj representation the den- 



sity of electrons is 



(36) 



A standard example of the theory which exhibits the 
topological instability and therefore superconductivity is 
the Dirac Hamiltonian 

H = i>l a { - iV + A) Vv + Ao^A/v - p^ a . (37) 

where a. = (a x ,a y ) and (3 are 2x2 Dirac matrices: 
{a x ,a y } = 0, (3 = —ia x a y . This Hamiltonian has 
5=- J F(r)dr zero modes (the flux is directed up). Wave 
functions of zero modes are: 



$( r ) =e -ifMr')-dr'-0fA(v') 

/3$o(z) = -$„(*) 



X dr 



$o(z), (38) 



where $o is any polynomial of degree ^- / F(r)dr - lE3. 
The fact of existence of zero modes implies that the en- 
ergy at p = has a linear term ( |l9| ) in chemical potential 
K(p = 0, r) = 2i5(r). Other models pf.rtepological super- 
conductivity were discussed in Refs.tali-l 

Let us comment on the relation between topological 
mechanism of superconductivity and superconductivity 
in the system of anyonsEilij. The models become very 
close after projection onto the low energy sector, where 
the relation (£0|) is treated as a constraint rather than 
as a result of minimization of the energy. The pro- 
jection can be done by introducing a Lagrangian mul- 
tiplier Aq for the relation (EG) and commutation rela- 
tions [A x (r),A y (r')] — 2niK (r — r'), i.e. by adding the 
Chern-Simons term with the kernel K to the Lagrangian. 
In anyon model the kernel K(r) is replaced by 2S(r) (even 



at p ^ 0), so that the relation between topological charge 
and the density becomes local p(r) = 2 



2tt 



This sim- 



plification results in a generation of transversal electric 
currents or "internal" magnetic field by light or by inho- 
mogeneous electric charge. Due to the same reason, the 
Meissner effect in the anyon model per se exists only at 
zero frequency, zero momentum, zero temperature, in- 
finitesimal magnetic field etc. These unphysical conse- 
quences originate from the topological constraint (p3). 
To avoid them one must determine the kernel K(p,r) 
self-consistently. 

The topological mechanism of superconductivity takes 
place in certain models of 2D doped Mott insulators. We 
point out several steps involved in the derivation of the 
Dirac Hamiltonian ( p7| ) in Appendix, while here we would 
like to give a physical interpretation of the vector poten- 
tial A and its flux F. The vector potential Ai is the 
phase of the hopping amplitude of an electron in direc- 
tion i = x,y. Its flux is the phase of the total hopping 
amplitude along a closed path. On a square lattice at 
small doping it is 



F(r) 



Im In trW, 



(39) 



where the chirality W is defined as 

W = (l+,r.S(r))(i+,T.S(r + e :c )) 

x (- + er ■ S(r + e y + e x )) 

and r,r + e x ,r + e x + e y are three consecutive points of 
a lattice and u are auxiliary Pauli matrices. 



(40) 



C. Hydrodynamics from the mean-field 
approximation 

The hydrodynamics of a superfluid ( |3l"| ) can easily be 
obtained from the model ([37]). Let us see how the en- 
ergy (|l9|) of a spin singlet state changes under smooth 
variations of "electric" E = dt A — 'VSp and "magnetic" 
F = V x A fields. To keep track of spin variations 
we add an external field A 3 to the Hamiltonian ( p7| ) 
zV — A — * iV — A — A 3 U3. In Gaussian approxima- 
tion the result consists of two separate parts — spin and 
charge. In the Coulomb gauge VA = the density of 
energy is 



SE = SE C + 5E S , 

6E c = ^ + ^F* + 5p(§-V «-,>»• 



SE, 



2 

n„ 



(41) 
(42) 



3\2 



(d t A 6 ) 



^-(VxA 3 ) 2 + ^A 3 xd t A 3 , (43) 



where polarization operators 

n = ^- 2 o-||(^,k)i||(-^,-k)>, 

= k~ 2 (j ± {Lj,k)j ± (-uj, -k) - j ll (uj,k)j ll (-uj, -k)), 
K = iu- 1 {j ll (u,k)j x (-Lj,-k) (44) 



G 



are transversal current-current correlators of free Dirac 
massive fermions at k, lo — > 0. Since the particles are 



massive the propagators (44) are known to be ITo = 
n_i_V(7 2 ~~ > const and K — > 2 . 

To obtain the hydrodynamics ( |30| ) one must rewrite 
this result in terms of spin and charge displacements. 
Minimization over <5/i gives the relation (BO) for the 
charge sector. Substituting this relation into ( |42] ) gives 
the hydrodynamics of the ideal liquid (|3l|). In the spin 
sector (^) the story is different. The total spin is kept 
to be zero. Therefore, the Chern-Simons term remains 
in the spin sector and gives the hydrodynamics of an in- 



compressible spin liquid, 
obtain the eq. (B2). 



Writing d t u s = —SC/SA, we 



D. Electron as a composite object. The vertex 
operator 

The most difficult part of the theory is to find a rela- 
tion between the true electron operator c a and "spinned" 
fcrmion tfj a . The problem is that the local field theory 
( |37j ) , is written in terms of a not gauge invariant operator 
Moreover, without gauge field (i.e. without interac- 
tion) ?/v being a Dirac spinor has 1/2 - orbital momen- 
tum. The electronic states are gauge invariant, and spin 
singlet states must have an integer orbital momentum. 
Another problem is that an electronic excitation carries 
a typical momentum of the order of fc/, while typical 
momenta of Dirac particles are close to zero. These dif- 
ficulties do not occur while studying the hydrodynamics, 
but rise in matrix elements. Below we employ two ap- 
proaches. In Appendix we derive the Dirac theory ( 1371 ) 
from a microscopical model of a doped Mott insulator 
and find the electronic operator and the vertex operator 
on a regular basis. In this section we conjecture the form 
of the electronic operator based on plausible physical ar- 
guments. 

The requirements for the electronic states are 

(i) electron is gauge invariant, i.e. remains unchanged un- 
der a non-singular gradient transformation A — > A+ VA; 

(ii) In a sector with completely filled zero modes, i.e. the 
flux and the number of particles obey the topological 
constraint (^), a charged singlet excitation is a spatial 
scalar, i.e. its wave function has a zero orbital moment 
I = 0; 

(iii) Since an electronic liquid is compressible, the 
most essential electronic modes have momenta k ~ 
(2^/p) 1/2 = k f . 

We find the electronic operator in three steps. First we 
take care of the gauge invariance. We define the vertex 
operator V a (r) as an operator which creates a flux quan- 
tum and a zercupode at the point r in the state with the 
spin a = ±l/2tL3. It obeys two relations 



V^ai-iV + A)V„ = -iaV, 
[F(r'),K(r)] = 2nV (T (r)S(r - r') 



(45) 
(46) 



These conditions are valid in the subspace of zero modes 
( p8|) and can be solved by virtue of the commutation 
relations of the Chern-Simons theory [A x (r), A y (r')] — 
i7T<5(r — r'). In terms of the density of particles with spin 
a: c\c a — ~ V • u a , or in terms of their displacements u CT 
the vertex operator can be written as 



VJr) — e 2wi u a {r')xdr' -2irf3 J* u„(r')-dr' 



(47) 



The operator Vj tp a is gauge invariant. 

The second step is the rotational invariance. Let us 
consider a wave function of the free Dirac field in two 
spatial dimensions. In the basis, where a- matrices are 
a x = a 3 , OLy = —<J2tP = &x, the solution of the Dirac 
equation with momentum p is e lpr e~i^ :iTgf - p '(up, v p ) 
for the positive energy E = +E p = ^Jp 2 + Aq, and 
e jpr e -f /3arg(p)^^ —Up) for E = —E p , where arg(p) is 
an angle of the momentum p, relative to the x-axis and 

u v = J h( 1 + B;) and v p = \f\( 1 ~ W> are the BCS 

wave functions. The spinor carries an angular momen- 
tum / = 1/2. We unwind the Dirac field by a chiral 
rotation 



Vv(p) 



,|/3arg(p) 



Vv(p)- 



(48) 



This singular transformation has a clear physical sense — 
it projects the spinor wave function onto a direction of 
the momentum p. Indeed, the chiral transformation 
in two spatial dimensions, where (3 = —i/2 [a x , a y ] is 
equivalent to a spatial rotation of the momentum p by 
the angle arg(p) which aligns the momentum p along 
the x-axis of the coordinate system. Now, without topo- 
logical gauge fluctuations, the transformed operator is a 
spatial scalar. 

The third step is to boost the fermion to the Fermi 
surface. In the chosen basis the upper and the lower 

components of the Dirac field %j) a = (^ip!j~\ ip^^j corre- 
spond to the states propagating forward and backward 
along the vector p. To construct an electronic operator 
we shift the momentum of the upper component by the 
Fermi vector kf = kf^ directed along the momentum 
p: p^k = kf + p and the momentum of the lower 
component by — kf : p — > k 2kf = — kf + p 



c a (k f + p) 
Co- (-k f + p) 



|/3arg(k; 



Vv(p). 



(49) 



Here we used c a to indicate that the gauge field has not 
been taken into account yet. 

Formally, the chiral rotation (^) may be also un- 
derstood in the following way. The wave function of 
the Dirac field depends explicitly on the choice of a- 
matrices, i.e. on the choice of holomorphic coordinates 
in a plane. The chiral transformation aligns holomorphic 
coordinates relative to each point of the "Fermi surface" , 
i.e. sets up the momentum dependent a-matrices 

Q k( = e f/3arg(k f ) Q , e -f,3arg(k f )^ 
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In terms of c the free Dirac Hamiltonian describes 
isotropic backward scattering 

H = J {&4(k)c(k) + A (c f (k)c(k - 2k f ) + /i.e.) } dk, 

(50) 

where — v {k ~ &/)< 

Assembling all pieces we obtain the relation between 
the Dirac fermion and the physical electron in the sector 
of zero modes 

V a {r)ca{v) ~ J e- ikr e-i axs ^^f = - 1 \la--kf)dk (51) 

where ipa~ 1 ' is the Dirac field, projected onto the zero 
mode sector ( ^^-ipa = ( ^ J ipj 7 f3 ~~ lS> ). 



IV. MATRIX ELEMENTS IN TOPOLOGICAL 
LIQUIDS 

A. Vertex operators 

To compute matrix elements we will need an opera- 
tor algebra for the vertex operator V a and two addi- 
tional vertex operators of the charge and spin sectors 
V = (V T ) 1/2 (Vi) 1 ' 2 and V s = (Vf) 1 / 2 ^)- 1 / 2 . Below we 
choose a holomorphic basis where (3 — 03 is diagonal. At 
(3 = — 1 (the flux up) the vertex operators depend on the 
holomorphic part of the displacement u a = u a ^ x + u a ^ v : 



Tr( \ iir f u(z)dz 

V[w) ~ e J v , 
s (w) ~ e J y ' , 



(52) 
(53) 



where it = (u-f + uj.)/2; u s = (itf — uj.)/2. Due to 
the commutation relations (|2^) vertex operators obey the 
operator algebra: 



V s (z)c r (z') ~ ( )V2 C ( 2 / )v;(2 , )j 



L 



y s (z) Ci (z')-(^)- 1/2 C i(^)^W, 



(54) 



^K(z') ~ (—-Y/\ a ( z >)V(z), 

V(z)V s (z')~V s (z')V(z), 
V{z)V(z') ~ V(z')V(z). 

Below we also use a two-particle vertex operator /j.(u, w) 
which creates a soliton and antisoliton of the spin dis- 
placement 



n(u,w) = V s (u)V s -Hw). 
The operator algebra gives 



(55) 



t t (u'-u^w'-w) 1 ' 2 

Pl(u ,w)cl(u)c\(w) = 

(u — w) 

x c|(u)c|(iu)/*(u', u/) (56) 
and if the points u, v! and w, w' coincide 



ci(u)ci(w)fi(u,w) 



(u — w) 



H(u,w)c^(u)c l (w). (57) 



The physical meaning of the vertex operators is 
straightforward. The vertex operator of the spin sector 
V s (z) creates a soliton of the spin displacement u s and 
removes spin up from the site z. Indeed due to (p9[), we 
have 



[V s (z),Vu s {z')] =V s {z)6(z - z'). 



(58) 



In contrast, the vertex operator of the charge channel 
V(z) creates a flux of the gauge field (the spin chirality), 
but commutes with displacements 



[F(w), V{z)] = 2irV{z)6(w - z), 
[u(w),V(z)]=0. 



(59) 



Two composite objects V s (z)c^{z) and V s (z)ci(z) are 
singlets but carry electric charge. 



B. Matrix elements 

An electronic operator c a does not create an elemen- 
tary excitation. The ground state with an extra particle 
is created by a composite operator which creates a flux 
and places a particle into the core of the flux. Creating 
a flux also generates spin waves. They are gapful. Let us 
first assume that the gap Ao is very large and we can ne- 
glect spin excitations. Then the vertex operator V of the 
charge sector creates a flux and a zero mode which gets 

occupied by a particle created by the operator tpa~ ■ 
Let us start from the two-particle state with zero mo- 
mentum: 



\N + 2) ~ £<J<jl J dvdv' 



(60) 



x V{v)V{v')^-^\v)^-- 1) \v')\N). 

A physical meaning of the composite operator which 
creates the ground state with an extra particle is more 
transparent in terms of gauge invariant electronic oper- 

In these terms the operators 



ators c a ~ V a 1 ipf 



i) 



(/?=-!) 



y s c T and V~V 



(/3=-l) 
J. 



V s 1 C| consist 



of the electron and the vortex (antivortex) of the spin 
density only. In its turn, the two-particle state (|60|) is 
composed of the two electrons and the operator Q55|), 
creating a vortex and an antivortex of spin density at 
the points of electron insertions. Electrons with opposite 
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spins see each other with vortices of the opposite angu- 
lar momenta I = ±1 (compare to Ref.EJ). This gives an 
angular momentum I = 2 for the pair. 

With the help of the operator algebra (^J,^,^) we 
are able to compute the tunneling amplitude as well as 
some other matrix elements. 

The singlet two- particle matrix element A(r) = 
(N\ Cl (r)c L {0)\N + 2) has the form 

A(r) ~ J dke- iarg(k V kr 

x{N + 2\V s (r)V s - 1 {0)\N + 2). 

Thus the two-particle matrix element is given by the 
correlation function of the vertex operators of the spin 
channel. The operator algebra V s (z)V~ 1 (z') = —3-7 : 
V s (z)V- l {z') : implies (N + 2| (r) V;- 1 (0) |iV + 2) ~ 

1 e - la rg(r) SQ that 
kfr 



A(r) i_ e -^rg(r) 

kfr 



e - l arg(k) e ikr dk 



(61) 



The matrix element ( |6l| ) is obtained in the limit of a 
very large gap Ao — > 00. This approximation is sufficient 
in order to analyze the transformation properties and 
the angular dependence of the tunneling amplitude. If 
the gap is not very large, an embedding of two extra 
particles creates a gapful spinon-antispinon excitations. 
These excitations do not interact with the zero mode of 
the charged sector and their wave functions (at kr 3> 1) 
are just the wave functions of an unperturbed theory (pC 



<&(+>(k, r) = u fe _ fe/ e* kr + Ufc _ fc/ e i < k - ak '> r , (62) 
$(")(k,r) = Wfe _ fe/ e lkr - Uk.k^- 2 ^. (63) 

Here the first (second) function corresponds to a posi- 
tive (negative) energy. The wave function of a spinon - 
antispinon pair with opposite momenta is 



A B cs(k,r) = e l T u k -k f Vk-k f 
= sign(/fc - kf) $ (+) (k, r)$(")*(k, r) 



(64) 



Ao 



vS + A 



Let us notice that in contrast to Cooper's pairing mech- 
anism, the gap in the topological mechanism of supercon- 
ductivity is generated via backward scattering and is of 
the insulator nature. Nevertheless, the wave function of 
a singlet spin excitation with zero relative momentum is 
the same as the BCS wave function 

To obtain matrix elements for a finite gap, we have 
to replace the factor e lkr in eqs. (|6l|) by the spinon- 
antispinon wave function ( |64[ ) 

A(r)~e- iarg(r) -A- / e - iarg(k) A BC s(k,r) ( ik. (65) 
kfr J 

A state with one extra particle is characterized by two 
momenta. One of them is a momentum k of the spin 



excitation and it is of the order of the Fermi momentum. 
Another one is a small momentum q of the zero mode. 
Arguments similar to the ones used in the derivation of 
(H) give 



(iV|c T (r)|iV + l,k,q) ~*«(k,r) 

/p -iarg(r')-iarg(k) 1 


V kfr' 



,iq(r-r ) i 



dr' 



(66) 



In momentum representation the one-particle matrix el- 
ement is 



(iV|c T (k')|iV + l,k,q)~ f dr'e^ r '( 



-i arg(r 7 )— % arg(k) 

kfr f 

x (u k _ kf 5(k' - k - q) + v k _ kf 6(k' - k q 2k,)) (67) 

These results lead to a number of important conse- 
quences. 



C. Orthogonality catastrophe and angle dependence 
of the tunneling 

(i) The angular dependence of the pair wave function 
and of the tunneling amplitude. The pair wave function 
consists of two vortices of the same charge — one comes 
from the wave function of the zero mode of the flux cre- 
ated for a spin up particle by a spin down particle (the 
factor e" larg ( r )f in @). Another vortex ( e " iar g( k )) is 
located at the center of the Fermi surface. The angular 
dependence of the pair wave function becomes obvious 
after the transformation: 



A(r) 



2arg(r)^w 



^iO ^ikr cos 9 



A,, 



kdkd6, (68) 



where 9 is the angle between k and r. 

The pair wave function forms the d-wave (/ = 2) irre- 
ducible complex representation of the rotational group. 
Similarly the one particle matrix element carries / = 1 
orbital moment. 

(ii) Orthogonality catastrophe. The overlap between 
two ground state wave functions with N and N + 2 elec- 
trons does not vanish in a large system. 

This is not the case for a one- (or any odd number) 
particle matrix element. An attempt to embed a sin- 
gle electron in the system leads to a half-occupied zero 
mode state. As a result this state turns out to be al- 
most orthogonal to all other states with the same spin 
and number of particles. The matrix element ( |67| ) van- 
ishes at q -> as (N\c a (k + q)\N + 1, k, q) ~ q 1 * 2 . This 
is a well-known phenomenon in electronic physics in one 
dimension^. A new feature is that it vanishes also as a 
result of the phase interference over the Fermi surface 
(the integral over argr in (|66|)). 







Due to the general arguments of the previous section, 
the tunneling is determined by the instantaneous pair 
matrix element (|65|). 

(hi) Tomographic representation. Eq. ( |68|) consists of 
the integral over the Fermi surface and may be viewed as 
a "tomographic" representation of the matrix element 



The result gives the universal dependence of the pair 



A(r) -e- l2arg(r ) <£dkD k (r) J A BCS {k,r)kdk 
where the propagator .Dk(r) is given by: 
Ac(r) 



r • k 



k 1 

_ _gj(argk-argr) 

lirxk r 



(69) 



(70) 



The tomographic representation in electronic liquids has 
been anticipated inE3c3. The new feature is that the 
relative phase of electron pairs at different points of the 
Fermi surface are correlated by the factor exp(zargk). 

(xv)Bremsstrahlung. It is instructive to rewrite eq. ( pq ) 
in momentum representation 



A(k) 



-il arg(k) 



A BC s(k-q)Ac(q)dq (71) 



where the propagator 
Ac(q) 



_ _J_ g i(argk-argq) 



q ■ k + iq x k q 

is a holomorphic function relative to k, and 

A 



A BC s(k) = 



V(6(k)-/^ + Ag' 



(72) 



(73) 



This representation clarifies the physics of the topological 
superconductor. An insertion of two particles in the spin 
singlet state with relative momentum k close to kf emits 
soft modes of transversal spin current with the propa- 
gator -Dk-k^q)- As a result of this: (i) ground states 
differing by an odd number of particles are orthogona£j; 
(ii) the BCS wave function is dressed by soft transver- 
sal spin modes. This is in line with ID physics and the 
bremsstrahlung effect of QED. The new feature is the 
phase of the matrix element of the soft mode ( fT2| ) which is 
the angle relative to the Fermi momentum of the electron. 
By contrast, in the BCS the density modulations are 
not individual excitations but are composed from Cooper 
pairs. The interaction between density modulations and 
pairs — the major effect of the topological mechanism — 
vanishes in the BCS. 

(iv) Momentum dependence of the two-particle wave 
function and tunneling amplitude. The momentum de- 
pendence of the amplitude of the pair wave function 
|A(k)| is drastically different from the BCS (§). In the 
vicinity of Fermi surface \k — kf \ <C kf the integrals j69| ) 
can be computed: 



lA^HA^I-^lsinh- 1 ^- 

£/ fc^ A 



(74) 



v(k-kf)/A v\k-k f \<.A 
sgn(fc - k f )\og(v\k - fe/|/Ao) v\k - fc/j > A 



wave function on a • ^he cons tant |A(fc/)| is not 
universal and depends on states far away from the Fermi 
surface. 

In contrast to the BCS gap function, the pair wave 
function (see Fig. |^) is asymmetric around the Fermi sur- 
face. It is peaked at scale 



A 



v(k - k f ) a e/ log 



which is much greater than A . This feature (already em- 
phasized in the Sec. [n] where we discussed the frequency 
dependence) is another consequence of the orthogonality 
catastrophe and is common in one-dimensional physicsc3. 
Also, the pair wave function has a logarithmic branch cut 
at v(k — kf) = ±iAo in contrast to square root singular- 
ity of the BCS function. This indicates that transversal 
spin current soft modes are emitted due to the tunneling. 



A(k) 





k r = 1000\ 








k r =10 



k/k f 

FIG. 2. The modulus of the pair wave function | A| is shown 
as a function of momentum k at Ao/wfc/ = 10 _1 , 10 -2 , 10 -3 . 
In contrast to the BCS gap function, the pair wave function 
is asymmetric around the Fermi surface (k/kf = 1). It is 
peaked at the scale v(k — kf) ^ e/ log -1 which is much 
bigger than Aq, and has a long tail away from the Fermi 
surface. 



V. THE PAIR WAVE FUNCTION IN THE 
COMMENSURATE CASE 

A. The model 

In this section we study the lattice model which can 
be derived from the canonical t-J model of doped Mott 
insulator under a set of physical assumptions. This set 
of assumptions and the derivation is outlined in Sec. VI. 

The model describes lattice fermions on a two- 
dimensional square lattice interacting with a gauge field. 
The number of fermions is close to the number of lattice 
sites and the electronic part of the Hamiltonian is the 
square of the hopping Hamiltonian 



h = m z 
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H = £ {|a)tf(a, b)(b| + \b)U(b, a)(a|} (75) 

(ab) 

with a fluctuating hopping amplitude. Here |a) is a site 
of a sublattice A of the square bipartite lattice, |b) is a 
site of the sublattice B, and (ab) denotes neighboring 
sites. The phase of a fluctuating hopping amplitude is a 
gauge field A 



U(a,b) = e iA * b . 



(76) 



The Hamiltonian ( |75[ ) describes the hopping of a particle 
within sublattice A and another particle within sublattice 
B 



B = l £ {V t (a)t/(a,a')V'(a') 

<aa'),<bb'> 

+ ^(b)C/(b,b'Mb')}. 



(77) 



The hopping amplitudes describe two consecutive hops 
which leave electrons on the same sublattice — here a, a' 
and b, b' are nearest sites of sublattices A and B 

U(a,a') = J2 U(a,b)U(b,a'), 

(ab)A(ba') 

U(b,b') = U{b,a)U(a,b'). (78) 

(ba)A(ab') 

The spin of the original electron has been identified with 
a sublattice, say, particles with spin up live only on the 
sublattice A whereas particles with spin down live on 
sublattice B (see the Appendix). 

It will be explained in the Appendix that the electronic 
energy achieves its minimum if the gauge, field ^4 a ,b forms 
a flux 7r per plaquette (flux hypothesis^xHl): 



a+e T ,a+e T 



a+e,+e„,a+e.. 



+e H ,a 



-1. (79) 



We choose this as a mean field value for the gauge 
field A a b and consider fluctuations around this mean 
field. The fluctuations are smooth and slow and can be 
treated in the continuum limit. 

In the mean field flux the hopping amplitude along 
the diagonal of the lattice cell vanishes U(r,r + e x + 
e y ) = due to the maximal interference of the hopping 
along two different paths connecting points r and r + 
£x + e y . Then the dispersion becomes e(k) ~ cos 2k x + 
cos2fca; and the Fermi energy e(k) = /i at small chemical 
potential consists of four small pockets centered around 
Dirac points k/ = k± t ± = (±^,±|) (see Fig. |l|). 

Let us therefore decompose the fermionic operator 
into four smooth movers: 



Vv(r) 



(80) 



In the following we refer to the smooth functions V'k/.o^r) 
as the continuum part and to the factors e* k/ r as the 
lattice part of the fermion operator Vv(r)- 



In this basis the mean field Hamiltonian can be written 
in the continuum limit as the square of the Dirac operator 



H — D 2 — (a x id x + a y id y ) 2 , 



(81) 



where 4x4 Dirac matrices {a x ,a y } — act on four 
fermionic modes ip^. = ip±±. The choice of these matri- 
ces (a gauge freedom) corresponds to a relabeling of the 
Dirac points and is limited by the symmetry group of the 
Fermi surface (there are only four different gauges). We 
choose Dirac matrices to be 



T"3 <X> T 3 , 



n <8> t 3 , 



(82) 



where r are Pauli matrices and the first (second) matrix 
in direct product acts on the first (second) index of the 
ijj±±. This choice of matrices corresponds to the lattice 
Landau gauge and the mean field lattice hopping ampli- 
tudes in this gauge are: 



a* 



U(r,r + e x ) = 1, 
U(r,r + e y ) = (-!)*. 



(83) 



The amplitudes of clockwise and counterclockwise diag- 
onal hoppings have opposite sign: 

i/3 — a x a y — [7(r,r + e x ,r + e x + e y ) (84) 
= -U(r,r + e y ,r + e x + e y ) = t 2 ® 1 = -(-l) x '. 

The gauge was chosen such that aj — —oci- The ad- 
vantage of this gauge is that electronic modes c CT] ±± are 
indeed slow — the mean field Hamiltonian (^lj) achieves 
the minimal energy at zero momentum. 

Fluctuations of the gauge fields around the mean field 
value can be easily incorporated into the mean field 
Hamiltonian (|8l|). Setting ip(a) — tp^ and ip(b) = ip±, 
we obtain: 



H 



1 

2m 



^{(iV-A) 2 +/3F}^. 



(85) 



Let us note that the diagonal hopping (the last term in 
(|35l)) exists only due to fluctuations. For a more detailed 
discussion see the Appendix. 

The Hamiltonian ( p5| ) together with the correspon- 
dence between continuum fields and the lattice ( p0| ) is 
the field theory for the doped Mott insulator on a bi- 
partite lattice. To make this theory complete one has to 
add the relation between fermionic fields ip a an d ori ginal 
electron c„. This will be done in section Sec. VB (sec 



(ED)-. 

Doping stabilizes the flux phase by causing the cre- 
ation of the flux of the gauge field in order to absorb 
extra electrons (see (|20|) ) . This, however, increases the 
magnetic energy due to chiral stiffness (unless ilie un- 
doped Mott insulator is a chiral antiferromagnctEj) . We 
assume that at a certain doping the magnetic energy loses 
the competition (if any) with the gain in electronic en- 
ergy. Then a doMdifi^x phase becomes a topological 
superconductorc§E 
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B. Lattice zero modes and lattice vertex operator 

To clarify the correspondence between continuum limit 
and the lattice, let us write the ground state wave func- 
tion of the mean fiel d Ham iltonian (81) explicitly. As 



discussed in the Sec. (HID), it is an eigenfunction of the 
matrix (3 (|84| ) with eigenvalue —1, if the flux is directed 
up (i.e., we are adding electrons). This state is doubly 
degenerate: 



zk i 



+ i(e i 



+ e 



f A = e lk ++ r 

Sb =e jk ++ r -e jk -- r + i(e lK - * -< 



ik_L _ r 



(86) 



Solutions are chosen, such that each of them exists on 
only one sublattice: /a = on sublattice B and /b = 
on sublattice A. The state f A is occupied by an elec- 
tron with spin up, whereas the state fs is occupied by 
an electron with spin down. These functions are not 
gaug e invariant (i.e. depend on the choice of the matrices 
( p2] , p3| , p4| ) ) . To construct a gauge invariant wave function 
of the electron one must multiply it by the lattice Dirac 
tail: 



* CT (r,C) = (7(r,C)/ CT (r), U(r,C) 



c 



(87) 



where the product goes over some lattice contour C . The 
wave function explicitly depends on the contour. 

Then the electron operator consists of the vertex oper- 
ator V a {r) (see (E^)), the spinned fermion operator and 
the lattice wave function 



c CT (r) =V- 1 (r)V a (r,C)Mr)- 



The factor ty a (r,C) may be considered as the lattice 
counterpart of the vertex operator. The form of the elec- 
tronic operator d87j , |S8| ) is to be compared with similar 
equations ( [38|]5l| ) — by decomposing into Fermi modes 
e ik ±± r we obtain a discrete analog of the phase factor 
e^/ 2 of ©. 

The theory (|85|-|s"8|) together with a correspondence be- 
tween continuum and lattice fields ( |S0| ) describes the elec- 
tronic topological fluid of the doped Mott insulator on a 
bipartite lattice. 

Below we will use the two-particle singlet zero mode 
state: 



vf(r,r') = U{v,C r )U{v',C' r )f A {v)f B (v'). 



(89) 



It depends on two strings (contours) ending in points 
r and r'. Fluctuations of the strings are physical ex- 
citations of the pair (not an artifact of the approach). 
The presence of string degrees of freedom is a general 
feature of gauge theories and becomes important when 
the coupling with matter is strong. In the commensu- 
rate case strings fall in four groups within which $ is 
the same. These groups correspond to the states with 
pairing from different Fermi points kf and ky, i.e. to a 
pair with a total momentum P = k/ + k* = (0,0) = 



(±7r, ±7r), (±7r, 0), (0, ±7r). The ground state pair wave 
function obviously has zero momentum P = 0. A class 
of strings which gives zero momentum to a pair is repre- 
sented by two contours following each other from some 
reference point up to a point r = (x, y) and then a single 
string along the y-axis to (x, y') and then to the point 
r' = (x',y') along the x-axis. In the chosen gauge this 
factor is: 



U{r,C r )U{r',C' r ) = (-lffo-V). 
Then the wave function is translation invariant 
*(r, r' 



sin ^ (X + Y) + i sin -{X - Y) 



E «- 

£/ =(±f.±f) 



iarg(k/) i k/R 



(90) 



(91) 



where arg(k/) = f , , ?f and R = r - r' = (X, Y). 

This is a direct discrete analog of the factor e tff +i- krcose 
in Eq. (|68|). The two-particle wave function forms an 
irreducible complex p-wave representation of the crys- 
tal group: under mr/2 rotation it produces the factor 
e -in7r/2 for n = 1, 2, 3, 4. It can be viewed as a discrete 
vortex located at the center of the Brillouin zone. 



C. The pair wave function 

Let us now proceed with matrix elements, 
culations are very similar to those of Sec. 



The cal- 



IV B, and we 



review them briefly. Let us add two particles in a sin- 
glet state with momenta kf ± p symmetric with respect 
to the center of the small Fermi surface or equivalently 
with momenta ±(k/ + p) symmetric with respect to the 
center of the Brillouin zone. The ground state with N + 2 
particles has the form: 



\N + 2)=e 



E 

=(±f,±f) 



ik,(r-iO f dvdv > dp 



*(r,r')A B cs(p,r-r') 
M (r,r')4(k / + p)4(-k / -p) \N), 



(92) 



where operator /i(r,r') — the creation operator of the 
spin vortices with opposite momenta — is given by (|5E|), 
\f(r, r') is the wave function of a singlet in the zero mode 
state (|9l]), and Abcs(p,i") is the BCS wave function 
d64|): The integral over p goes over a small Fermi sur- 
face (pocket). Then the pair wave function in terms of 
Fermi components (|80| ) has the form: 

A(n - r 2 ) ~ (Ar|c T (n) Ci (r 2 ) f dvdv ' 

fe/=c±f ,±§r 

*(r,r') ABcs(k/,r-r'V(r,r')4(k)4,(-k) \N) 

giarg (n— r 2 ) 



*(ri,r 2 ) 



E 



ikjr(ri-r 2 ) 



k /= (±f ,±f) 

x AbcsOc/.ti -r 2 ){N\(i(n,r 2 )\N). 
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Using eq. ( |9l| ) we obtain the pair wave function in two 
equivalent forms 

, sin? (X + Y) +isinf (X -Y) , s , s 
A(R) ^ ^ + iY 21 ^Abcs(R) (93) 

e- larg(k ^e ik ^ R A B cs(R), (94) 



X + iY 

fc /= (±f,±f) 

where R = r — r' = (X, Y). The numerator of the first 
expression (^3|) is the discrete analog of the continuous 
holomorphic function in the denominator. Under -k 12 ro- 
tation it produces the factor e~~ l7T / 2 . Another e" l7r ' 2 fac- 
tor is produced by the continuum part. Both phases add 
to e _i7r = —1, so that the tunneling amplitude belongs 
to the irreducible d-wave representation of the group of 
rotations of square lattice 



A(-Y,X) 



-A(X, Y). 



(95) 



The second expression ( |94| ) is the commensurate ver- 
sion of the cq. (pq). Let us look at the tunneling ampli- 
tude in momentum space. It is: 



A(k) = £ 

fc/=(±f ,±\ 



3 - i arg(k / )- l arg(k-k / )j(| k _ (gg) 



where arg(kf) = f, 4p X' T an< ^ f(p) * s a smooth func- 
tion. The tunneling amplitude consists of two similarly 
oriented vortices: one in the center of the Brillouin zone 
(lattice part) while another is at a Fermi point. Both of 
them contribute to the phase (Fig. pi) 



A( ky, k x ) — A(k x , ky). 



(97) 



The pair wave function in both commensurate and in- 
commensurate (eas. (|7l|) , (p6|) ) case may be written in a 
unified form (|5|jg,0) with a g enera l dispersion e^. In the 



incommensurate case of Sec. IV C| e(k) — /i ~ «o(|k| — kf) 



and the Fermi surface is a circle. In the commensurate 
case of the Sec. M e(k + (±tt,0)) = e(k + (0,±tt)) and 
has a maximum (minimum) at k = (±tt/2, ±7t/2). If the 
chemical potential is close to the maximum (minimum) 
of e(k), the Fermi surface consists of four pockets with 
dispersion e(k) - n ~ Ek /= (±f ,±f ) u o(l k | - kf). 



VI. DISCUSSION 

1. Electrons dressed by soft modes of density modu- 
lation. The models of topological fluids have a lot in 
common with electronic physics of one spatial dimen- 
sion. At the core of this physics are zero modes (or ax- 
ial current anomaly) and the orthogonality catastrophe. 
Neither electrons nor Cooper pairs are elementary exci- 
tations in topological fluids. The one-electron (any non- 
singlet state) insertion drastically changes the ground 
state of the system, so that the matrix element between 
two ground states with N and N+l electrons vanishes in 



a macroscopic system. The matrix elements of two parti- 
cles in a singlet state do not vanish but are significantly 
modified by the interaction with soft modes of density 
modulations. In particular, the poles of Green functions 
on the mass shell are replaced by branch cuts. All of 
this is known in one-dimensional physics, but now we see 
similar features in dimensions greater than one. One of 
the reasons for this is that the emission of a soft mode is 
asymptotically forward: it does not push an electron to 
another Fermi point (one can see that the phase factor in 
eq. (|7l|) suppresses the scattering channels different from 
forward scattering). 

Eq. (|7l]) clarifies the physics of a topological super- 
conductor. A pair of electrons emits soft transversal 
spin modes and intepwcts with them (compare with one- 
dimensional systemsc3). These soft modes also exist in 
the BCS theory, but their interaction with electrons van- 
ishes at the mean field level and is negligible beyond the 
mean field. 

2. Excitations are solitons in charge and spin sectors. 
Similar to one-dimensional electronic liquids the elemen- 
tary excitations of topological electronic liquids are soli- 
tons of the spin and charge distortions u s and u. They 
carry spin 1 /2 and no electric charge and charge 1 and no 
spin respec tively. The vertex operators V s and V of the 
Sec. ( IV A ) are the creation operators of these particles. 
An interaction between solitons of different sectors van- 
ishes at small momenta, so the spin and charge channels 
are practically separated. 

3. Complex d-wave pair wave function with no nodes 
in the gap function. A new feature of two spatial di- 
mensions (in addition to the very existence of supercon- 
ductivity) is that not only the amplitude but also the 
phase of the two-particle matrix element is strongly mod- 
ified. The pair wave function acquires a phase which is 
given by A(r) ~ e~ 2largr |A(r)| and realizes the complex 
d-wave representation of the group of spatial rotations. 
The physics of this d-wave state is drastically different 
from the conventional d-wave. The amplitude of the pair 
wave function as well as the gap function has no nodes. 

4. Asymmetric pair. The interaction with soft modes 
gives rise to a broad structure of the amplitude of the pair 
wave function in momentum space around what one may 
call Fermi surface. The width of this structure is of the 
order of fc//log(e//Ao) which is much bigger than Aq/v, 
the width of the peak of the BCS wave function. It is 
rather straightforward to take into account the frequency 
dependence of the two-particle matrix element. In order 
to do it one must replace Abcs and -Dk(q) in the eq. (171 
by 



A B cs(^,k) 



Ar 



£>k(w,q) 



lo 2 — v 2 q 2 



i(arg k— arg q) 



(98) 
(99) 



so that 
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A(fi,k) ~ e - 42ar s( k ) J dqdu 

x A BC s(« - w, k - q)£> k (w, q). (100) 

The orthogonality catastrophe and a strong interac- 
tion with soft modes (in both spin and charge fluids) 
also affect other matrix elements. In particular, the one- 
particle correlation function acquires a branch cut on the 
mass shell u> — yj (e(k) — /i) 2 + A 2 , rather than a pole. It 
leads to a broad structure around the Fermi surface in 
the photoemission spectrum, and to the distribution of 
the numbot, of particles, similar to those known in one 
dimensioned. It also affects the I- V characteristics of di- 
rect tunneling (0,||). 

The vertex operator algebra developed in this paper 
allows one to compute the long time, large distance be- 
havior of the most interesting matrix elements and cor- 
relation functions. 

5. Parity and time reversal symmetry breaking. There 
have been misconceptions in the literature regarding 
parity symmetry breaking of the ground state of two- 
dimensional topological liquid. The following comments 
seem to be in order. The spatial parity and time reversal 
symmetry are simultaneously broken in the ground state. 
This reflects the chiral nature of zero modes. However 
it is not easy to detect this symmetry breaking experi- 
mentally. The reason is that, although the time reversal 
symmetry is broken, there are no spontaneous local elec- 
tric currents neither in the bulk nor on the edge in any 
steady state. One can see it from the hydrodynamics of 
the charge sector (|3l]), but this fact remains valid even 
beyond low frequency range. Even more general, all diag- 
onal singlet matrix elements are parity and time reversal 
even. Therefore, one should not expect to observe the 
time reversal symmetry breaking by-measuring Faraday 
rotatiorcil and muon spin relaxations. 

Another matter are off-diagonal or non-singlet matrix 
elements. The broken time reversal symmetry is explicit 
in the complex d-wavc tunneling amplitude. Another 
manifestations of the broken symmetry can be seen in 
the spin sector. Among them is an expectation value of 
the spin chirality S ■ VS x VS and a novel feature — edge 
spin current. This is a two-dimensional version of the 
known phenomenon in ID. A spin chain with gapful bulk 
spin excitations develops gapless spin excitations at the 
edges. In two-dimensional spin liquid edge excitations 
are chiral spin currents. Edg£-piagnetic excitations have 
been observed in spin chainsEj. One may expect to find 
these soft edge spin excitations in model systems with 
an enhanced boundary (say an array of superconducting 
islands). 

Another obstacle to observe time reversal symmetry 
breakiufris that parity alternates between odd and even 
layersESEj. The genuine parity breaking takes place only 
in systems with odd number of layers. Any realistic junc- 
tion, however averages over many layers. The tunneling 
amplitude averaged over two layers is A ~ cos(2 arg k+a) 
where 2a is a relative phase between the layers. This 



phase is determined by the intra-layer crystal anisotropy 
and most likely is locked to be zero. In this case the phase 
of the pair wave function is indistinguishable from the 
conventional d-wave. In a hypothetical mono-layer tri- 
crystal type of experiment one would expect the trapped 
flux to be an integet, (in contrast to a half integer for 
conventional d-wave£5). 

6. In addition to the search for different manifesta- 
tions of broken time reversal symmetry one may try to 
test somewhat more modest predictions of the theory in 
cuprates. Among them are: (i) Magnetic excitations on 
the edge. Although the current averaged over layers van- 
ishes, one may try to detect softened magnetic modes 
on the edge; (ii) Asymmetry of the shape of pair wave 
function A(fc) and its broad structure (see Fig. ||). (iii) 
Absence of nodes in the gap function. To the best of our 
knowledge the upper bound foc-the minimal gap (in (110) 
direction) set by experimentst£l is about 3% of the gap 
in the crystal axis direction. This is consistent with the 
gapless spectrum but also leaves some room for a small 
gap. 
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APPENDIX: DOPED MOTT INSULATOR AS A 
TOPOLOGICAL FLUID 



The phenomenological picture described in Sec. Ill 
may be realized in certain models of the doped Mott insu- 
lator. We start from the canonical Hubbard model with 
an infinite on-site Coulomb repulsion or the t-J model on 
a square two-dimensional lattice 



H 



where 



% C ia- C J CT Jij^i 



(Al) 



(A2) 



is the spin operator of an electron at the lattice site i. 
The hopping amplitude and the antiferromagnetic ex- 
change amplitude Jij > connect the nearest neighbors, 
the total number of electrons is close to the number of 
lattice sites: N e — Nq{1 — 5), while a strong Coulomb 
interaction does not allow doubly occupied states: 
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E 



C ia C ia 



(A3) 



This model does not have distinct scales, which are nec- 
essary to separate the physics of topological fluids from 
a variety of other physical phenomena which occur in a 
correlated electronic system. We, therefore, are forced to 
proceed through the adiabatic approximation, which is 
not justified by any small parameter, but correctly cap- 
tures the physics of interest. In fact the simple model 
( Al) can be generalized to make the approximation para- 
metrically correct. We will not concentrate on this, but 
rather carefully emphasize a set of assumptions. In the 
adiabatic approximation one considers the hole motion 
in a slowly varying spin configuration. In this case we 
may use a semi-classical strategy: first find a static spin 
configuration which minimizes the energy of the system 
with a given doping, and then take into account quantum 
fluctuations around the static mean field. 



1. Adiabatic approximation 

a. Hopping amplitude, chirality operator 

Adiabatic arguments may be applied to a model where 
fast processes are integrated out. In the t-J model as it 
is, the dynamics of spins is not adiabatic. The reason 
is the short distance antifcrromagnetic spin correlation. 
Each jump of a hole abruptly changes the spin config- 
uration by flipping a spin on a sublattice and increases 
the energy of the system by A,/ ~ J. However, two con- 
secutive jumps bring a hole to the same sublattice, so 
that a spin exchange between two shifted spins may heal 
a wounded antiferromagnetic bond. As a result the spin 
configuration remains approximately unchanged only in 
the second (even) order in the hopping process. 

To estimate an effective hopping Hamiltonian, in which 
single jumps (odd number of sites) are integrated out (a 
sort of Schrieffer- Wolff transformation), we assume that 
the kinetic energy of the hole is smaller than the exchange 
energy of the wounded antiferromagnetic bond (the anti- 
ferromagnetic correlation energy). It allows one to con- 
sider the hopping from, say sublattice A to sublattice B 
as a virtual process. In other words t he m atrix elements 
of the hopping part of the t-J model (Al) at low energy 
states are small — an effective Hamiltonian occurs as a re- 
sult of the second order of perturbation theory in t. A 
quantum spin liquid with a spin gap Aj, where spin- flip 
excitations are separated from the ground state, provides 
favorable conditions for the adiabatic approximation. It 
is required that tS <C Aj for the validity of the adiabatic 
approximation. 

Let us consider a matrix element of the hopping am- 
plitude between nearest sites a and a' of the sublattice 
A in the second order of t. Let \aa) be the low energy 
states of the antiferromagnet with site a removed and 



with a being the spin of the "removed" electron. Then 
the hopping amplitude is: 



A aa ,(a,a') =t 2 ^^] 

aa' \b) 



{aa 



\b)(b\c b0 cl, \a'a') 



(A4) 



where the sum goes over intermediate states \b) with site 
b removed and Eo (Eb) is the energy of the ground (ex- 
cited) state. 

The following three approximations follow from strong 
short range antiferromagnetism. 

(i) First of all one can replace \aa) by c acr |{S}), where 
|{S}) is a low energy state of the undoped antiferromag- 
net: the spin configuration of the ground state of the 
antiferromagnet with and without removed site is not 
drastically different in the vicinity of the removed site. 

(ii) Next, we may replace Eq — Eb by — Aj and take 
the denominator out of the sum in (A4). Indee d, i nter- 
mediate states \b) which contribute to the sum (A4) are 
different from the ground state by a permutation of spins 
on sites a and a' . Their typical energy is of the order of 
A j. Matrix elements of other states with low energy 
vanish. This means that the time which a hole spends 
on the sublattice B is very short, so the two consecutive 
hopping operators in (A4) act at the same time. As a 
result, 



<{S}| \b)(b\ <|{S}> 



\h) 



Eq — Eb 



t 

~ -^-Eh <{ S }l C L C a« C L C b/3 C l'/3 C a' CT '|{S}), (A5) 

where the sum in the last equation goes only over sites 
that belong to a two step contour connecting points a 
and a' . 

The effective hopping amplitude may be re-expressed 
entirely through-spin operators, namely through the chi- 
rality operatort£l 



t 



A aa , (a, a') = — ^ Woo> (a, b, a'), 



(A6) 



where W(a, b, a') is the chirality operator which drives a 
hole around a closed path a — > b — > a' — > a 

W(a,b,a') ee ({S}\ct a c a J ba c b/3 cl c a , a ,\{S}) 

= (i +< r.S o )(i+«r-S t )(i+<r.S 8 0. (A7) 



The effective Hamiltonian, then has the forr 
H= {4(a)A^(a,a')c.'(a') 

(aa').(bb') 

+ 4(b)A CTCT ,(b,b'K,(b')} 
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+ A./{4(a)<x CTCT ,£v(a) ■ S 6 + 4(b)o- CTCT /c CT /(b) • S Q } 

(ab> 

+ J s °> ■ S b , (A8) 

(ab) 

where a, a' and b, b' are the nearest points of sublattices 
A and B respectively. 

b. Gauge fields 

The hopping Hamiltonian ( |A8|) is not of practical use — 
the hopping amplitudes W aa i (a, b, a') and electronic op- 
erators are not independent. The practical way is to in- 
troduce a gauge field and a non gauge invariant fermionic 
operator instead of the electronic operator. Then the 
phase of the hopping amplitude will be a flux of the gauge 
field. Let us rotate all spins to the third axis by the U (2) 
matrix g 



a ■ S = g^a 3 g 

and introduce fermionic operators 

"0(a) = ST<xC CT (a), 
^( b ) = AW (a) 



(A9) 



(A10) 



on sublattices A and B. Then the hopping Hamiltonian is 
expressed through a non gauge invariant but independent 
fields ip and U(a, b) 



(AH) 



H = Yl ^ f (a)C/(a, b)V(b) + h.c. 

<ab) 

where the flux of the field U(&, b) is the chirality 
W(a, b, a') = trW = [/(a, h)U(b, a')U(a', b) 

= i + is Q • S b + is a • S„, + is Q , ■ S h (A12) 

+ iS a x S 6 ■ S a / 

This is the effective lattice Hamiltonian of the doped 
Mott insulator. 



c. Continuum limit: The flux phase stabilized by doping 

The next step in the adiabatic scheme is to separate 
the fast fields from the slow ones and then proceed to 
the continuum limit. To do this let us first d eter mine 
the mean field value of the hopping amplitude ( |A6| ), i.e. 
a static spin configuration which minimizes the the en- 
ergy. Among a variety of local minima we concentrate on 
a liquid state, and discard charge (spin) density waves 
which break the crystal symmetry. The most interest- 
ing picture appears if the electronic part rather than the 
exchange part of the Hamiltonian dictates the charac- 
ter of the magnetic state. This happens if the gain in 



renormalized kinetic energy of dopants is larger than the 
difference between the energy of the undoped antiferro- 
magnet (Neel state) and the true ground state. We as- 
sume that the true ground state of the doped antiferro- 
magnet belongs to the so-called flux phase universality 
class. This means that the kinetic energy of holes is es- 
sentially zero which gives a gain of the order of S 2 t 2 /J 
when compared with the Neel state. The difference be- 
tween the magnetic energy of flux and the Neel state Aj 
is of the same order but numerically less then J. Com- 
paring these energy scales we obtain that the flux state 
becomes energeticall y favor able starting at some critical 
doping S c ~ (J/t)^Aj/J. Although J > t the crit- 
ical doping may be— numerically small. Similar conclu- 
sions were drawn ii£3 although on the basis of a different 
model. This condition seems compatible with the adia- 
batic approximation. The latter requires that the typical 
renormalized kinetic energy of holes in the flux phase 
~ 5t 2 1 A j be bigger than the typical energy of spin exci- 
tations (spin gap). 

The moduli of hopping amplitudes A(a, a') and A(6, b') 
are determined by the competition between electronic 
and magnetic parts of the Hamiltonian. Due to the crys- 
tal symmetry: |A(r,r + e x )| = |A(r. r + e y ) | , and in 
the liquid phase they do not depend on the lattice site. 
We treat the mean field value of the modulus as a phe- 
nomenological constant and later set it to 1. However the 
phase of the hopping amplitudes may be inhomogeneous 
even in a liquid. The character of electronic processes is 
very sensitive to the phase. If the kinetic energy of elec- 
trons is greater than the magnetic energy, the phase must 
be chosen to minimize the electronic energy at a given 
values of moduli of hopping amplitudes. In other words 
the exchange part of the Hamiltonian directly governs 
the dot product of spins on different sublattices, while 
the cross product of spins is determined by electronic 
processes. .— ■ 

The flux hypothesis 2 -! suggests that the electronic en- 
ergy achieves its minimum if the chiralities along con- 
tours W(a, a + e x , a + 2e x ) and W(a, a + e y , a + 2a,,) 
are positive, while the relative phase between amplitudes 
along two different paths connecting sites on diagonals of 
a crystal cell 

W(a,a + e x ,a + e x +e y ) = \W\e^' 2 , 

W(a, a + e y ,a + e x + e y ) = \W\e-^' 2 , (A13) 

is the same in every crystal cell and is equal: 



$ = tt(1 - 5). 



(A14) 



The phase $/2 (i.e. the interference of different paths) 
is the flux through a clockwise oriented closed triangular 
paths within a plaquette. 



Let us note that the mean- field flux (A14) does not 
vanish as the doping S — * 0. We, therefore, set it equal to 
<f> = 7r and take the correction into account in the contin- 
uum limit. Then the hopping amplitudes along diagonals 



1G 



A(r, r + e x + e y ) — 2\W\ cos $ vanish in the mean field 
approximation. In other words, the two contributions to 
the diagonal hopping amplitude along two different paths 
have opposite phases and cancel each other. It means 
that chirality of two differently oriented triangular con- 
tours have opposite sign (W(a., a + e x , a + e x + e y )) = 
— (W(a, a + e y ,a + e x + e y )). The flux it per plaque- 
tte implies anticommutativity of translations of a holon 
r — > r + e x — > r + e T and r — > r + — >r + e x : 

U(t,t + e x )U(r,r + e y ) + U(r, r + e y )U(r, r + e x ) = 0. 

(A15) 

The Fermi-surface of the mean field state consists of 
four pockets around Dirac points k/ = k±,± = (±|-,±^), 
so we decompose the electron operator into four smooth 
movers 



Finally all fields are smooth and slow and we may take 
the continuum limit 



(A16) 



In the following we refer to the smooth functions c CTj ±±(r) 
as the continuum part and to the factors e lk/r as the 
lattice part of the fermion operator c CT (r). 

In this basis the mean field Hamiltonian can be written 
in the continuum limit as a square of Dirac operator 



H = tD 2 = t(a x id x + a y id y ) 2 , 



(A17) 



where 4x4 Dirac matrices {a x , a y } = act in the space 
labeled by the four Fermi points (±±). 



d. Fluctuations and the field theory 

Now we are ready to take into account smooth fluctua- 
tions of the phase of the hopping amplitudes (fluctuations 
of moduli are not that important): 



U(t, r + e<) = U(r, r + e»)e 



and proceed to the continuum limit. Then the hopping 
amplitudes along the crystal axes become: 

U(r,v + ei)U(r + e 4 ,r + 2e 4 ) = e ^*W+w«(>^) ) (A18) 

whereas diagonal hopping is only due to chirality fluctu- 
ations 

U(r, r + e x )U(r + e x ,r + e x + e y ) 
+U(r, r + e y )U(r + e y , r + e x + e y ) ~ F. (A19) 

The hopping Hamiltonian has the form 

H = -t{ ^(r + 2e 4 )e l ^ (r)+ ^ (r+e ' ) V'(r) 

i— x, y 

+ ^ (r + e x + e y )F(r)ij;(r) 

+ V f (r + e x - e y )F(r)i/j a (r) + h.c.}. (A20) 



H=^{{iV -A) 2 
2m 



(3F}ip. 



(A21) 



We find the hopping Hamiltonian to be Pauli operator. 
In order to compute matrix elements of the physical elec- 
tronic operators, one must complement the Hamiltonian 
by the relatio n between ip in the continuum and the c a 
on the lattice ( |A16| , |A10| , |A9| ). 
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